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Abstract 

We perform a Inonii-Wigner contraction on Gaudin models, showing how the integra- 
bility property is preserved by this algebraic procedure. Starting from Gaudin models we 
obtain new integrable chains, that we call Lagrange chains, associated to the same linear 
r-matrix structure. We give a general construction involving rational, trigonometric and 
elliptic solutions of the classical Yang-Baxter equation. Two particular examples are explic- 
itly considered: the rational Lagrange chain and the trigonometric one. In both cases local 
variables of the models are the generators of the direct sum of N e(3) interacting tops. 
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1 Introduction 



In the present paper we are concerned with an extension of Gaudin model. Since it was intro- 
duced in the 70's the Gaudin model jUIH] attracted considerable interest among theoretical and 
mathematical physicists, and played a distinguished role in the realm of integrable systems. Its 
peculariar properties, holding both at the classical and at the quantum level, are deeply con- 
nected with the long-range nature of the interaction described by its commuting Hamiltonians, 
which in fact yields a typical "mean field" dynamic. 

The algebraic counterpart of these physical features lies in the linear r-matrix structure, that 
characterizes both the classical and the quantum model, and holds whatever be the dependence 
(rational, trigonometric, elliptic) upon the spectral parameter. The linear r-matrix structure 
does not depend on the site number N (i.e. on the number of copies of the underlying Lie alge- 
bra) and moreover is not affected by a Inonii-Wigner contraction procedure. Such robustness 
property can for instance unable one to inherit the integrability property of the Lagrange top 
model from that of a two-site rational Gaudin magnet fTQj • More generally, it allows one to con- 
struct, out of the iV-site Gaudin model, novel interacting systems (both classical and quantum), 
whose underlying algebraic structure has been denoted in the literature as "jet extensions" of 
Lie algebras E] . 

In this paper we will elucidate the connection between Lie algebra contractions and Lie 
algebra jet extensions in a general and systematic fashion, and we will derive the associated 
families of classical (and quantum) integrable models. 

We have to mention here that pioneering results in this direction were obtained already in 
the 80's by Faddeev and Reshetikhin in their seminal paper ^2]) also, we have to mention the 
results obtained by Harnad in [J] for the rational case. However, we want to stress the connection 
among the systems investigated in the present paper and the Gaudin models, that to the best 
of our knowledge, is established here for the first time. 

2 Drinfel'd— Belavin solutions to cYBE and Gaudin Models 

Let V be a Poisson manifold whose symplectic leaves have generically dimension N. Differ- 
entiable functions on V form a Poisson algebra and a completely integrable system on V is 
specified giving an abelian Poisson subalgebra of dimension N. In the modern theory of clas- 
sical integrable systems the generators of this abelian Poisson subalgebra {Hi, . . . ,Hn}, called 
Hamiltonians of the system, are usually obtained as spectral invariants of a Lax matrix L(X) 
depending on a complex parameter A and whose matrix elements are functions on the Poisson 
manifold V. Let us denote with M the dimension of the square matrix L(X) defined by 

L{\) = Y,L ij {\)E ij , 

ij 

where (Eij)ki = 5ikdji is the canonical basis of the M x M matrices. Introducing 1 as the 
identity matrix of dimension M we define {X(A) ® 1, 1 ® L(fi)} as the matrix of Poisson brackets 
between the entries of L(X): 

{L(X) <g> 1, 1 ® L(jj)} = ^{LijiX), LMjEij ® E kl . 

ij,kl 

In this context a theorem due to Babelon and Viallet [2] ensures that the commutativity of the 
Hamiltonians {Hi, . . . , -ffjv} is equivalent to the existence of a M 2 x M 2 matrix ri 2 (X, (j,) such 
that: 

{L(A) ® 1, 1 <g> L(jj)} = [n 2 (A, n),L{X) ® 1] - [r 2 i(fx, A), 1 <g> L(jj.)], (1) 
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where 

»"2i(A,/i) = ilri 2 (A,/i)Il, 
with II the permutation operator on C M (g) C A/ : 

n(x <g> y) = y ® x. 

In this paper we will deal with an interesting subcase for equation (^Q) obtained when one 
imposes the following constraints on ri 2 and r-ix'- 

ri 2 (A,/i) = n 2 (A - /i) = -r 2 i(A - /i), 

and moreover requires that the r— matrix be "not dynamical" (i.e. independent of phase space 
variables). 

In this case equation Q becomes (omitting the indices 12): 

{L(A) (81,1® L{ji)} + [r(A - fi),L(X) <g> 1 + 1 ® L(/x)] = 0. (2) 

It is well known that, in the case of a constant r-matrix, a sufficient condition for the r.h.s. 
of to give rise to a well defined Poisson bracket is that r satisfies the so-called "classical 
Yang-Baxter equation" (cYBE): 

[ri 3 (A), r 23 (M)] + h 2 (A - m), ri 3 (A) + r 23 (At)] =0, (3) 

where 

n 2 = r (8)1, r 23 = l®r, n 3 = (1 ® IS) r X2 (1 ® II). 

This equation was studied by Drinfel'd and Belavin [Hj in the case when the r— matrix is 
an element of g ® g, with g being a finite-dimensional simple Lie algebra. Once one has such a 
solution r(A) of equation @, then (pi <g> p 2 )r(A) with pi and p 2 faithful representations for g, 
gives a matrix solution for the cYBE. 

Let {^"}*™ 1 S be a basis of the Lie algebra g, which is defined by the commutation relations 

x a ,x?] = cfx^ 1 , 

and let g a p be the inverse of the Cartan-Killing metric related to this basis. Drinfel'd proved 
that under suitable non-degeneracy conditions on r(A) the only solutions to cYBE © are given 
by: 

g aP X a ®XP 



rational solution: r(A) 



A 

oo 



trigonometric solution: r(A) = — — \ — . 

^ A — nuj 



n=— oo 



. elliptic solution: r(X) - V V"*- VDitoX- *X» 

^— ' A — nwi — rao)2 

n,m=— oo 

Here A and -B denote two finite order commuting automorphisms of g not having a common 
fixed vector. Since sl(n) is the only simple Lie algebra possessing such two automorphisms it 
follows that the elliptic solution can be defined only in the case g = $l(n). 

For our purposes, the relevant property of Drinfel'd solutions is that all of them can be 
generically written in the form: 

r(A)=^r®lY(A), (4) 
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for suitable choices of the functions f a (X). 

From the cYBE we obtain a system of functional equations that must be satisfied by the 
functions f a (X). In fact the first member of (jHJ) is given by 

'g a/3 X a ®l®X? f a (X),g lS 1®X^®X S / 7 (A) = 

= g aP g-ys Cf x a ® Xi® x^ f a (X) p{fj), (5) 

while from the second one we obtain 

g af3 X a ®xP®l / Q (A - /x), g iv (X 7 ® 1 ® X" / 7 (A) + / 7 (^)) 
= <M f a {X - M) [CT X" / 7 (A) + Cf X a ® X 5 ® X" / 7 (//) 

= <M 5 7 «5 ^ a ® X 7 ® X" Cf [/ 5 (A - M )/^(A) - f (A - , (6) 

where the last member is obtained from the previous one using the identity g^ v C" 7 = —g^s C^ 7 
and swapping the indexes 7 <-> /3, a «-> (5 in the first term and 7 <-> (5 in the second one. 
Equating Q and © we obtain the system of functional equations: 



dim g 



£ <fo/3 5 7 5 Cf [/ Q (A)/ 7 (/i) + / 5 (A - M )/*(A) - / a (A - p)f s (p) 

0,8=1 



(7) 



for all a, 7, 77 = 1, ...,dimg. Note that in equation (J7| the summation is not on all the repeated 
indices, so that the sum is indicated explicitly. 

To each matrix r(A) we can associate a Lax matrix whose entries are functions on the Lie- 
Poisson manifold V s associated with g in the following way: we choose a faithful representation 
p of g and we denote by {y a }^=i the coordinate functions related to the basis dual to {^ a }^™ 9 : 



} = Cfy\ 



(8) 



x a ,x p =cfx~< => 

We define our Lax matrix as: 

L(X) = T + g af3 p(X a )yl 3 r(X), 

where r E p(g) is a constant matrix chosen in such a way that 

[t® l + l®r,r(A)] =0. 

This Lax matrix satisfies the linear r-matrix algebra 

{L(A) ® 1, 1 ® L(p)} + [r(A - p),L(X) ® 1 + 1 ® L(p)} = 0, 

with respect to the ordinary Lie-Poisson bracket, where r is given by Q). 

Now let us consider the Lie-Poisson manifold associated with the direct sum of N copies of 
a simple Lie algebra g. We denote with {yf}^=i S the coordinates on the dual of the i— th copy 
of g: 

{y?,yl}=5 i0 Cfy]. (9) 
On each site we can define a matrix Lj(A): 



L i (X)=g a pp(X a )yPf a (X-X l ), 



(10) 
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where for any i we have shifted the spectral parameter by a fixed quantity A» € C. Each matrix 
Li will obviously satisfy equation © and thanks to the locality of the Poisson relations © so 
will do also the total Lax matrix: 

N 

Lg(X)=r + Y J U{X). (11) 

i=l 

It follows that the spectral invariants of the Lax matrix Lg(X) will provide us with a set of 
involutive Hamiltonians that define the N— sites Gaudin model associated to g [4"| IK} IT3"1 |Q| ITT]. 

3 Algebraic contractions of Gaudin models 

In the present section we illustrate our main result. Starting from Gaudin Lax operators 1)11(1 
we obtain, through a Inonii-Wigner contraction, new Lax matrices satisfying the same r-matrix 
structure pertaining to the starting ones. 

Let us recall the notion of generalized Inonii-Wigner contraction |Sj. Let V a complex finite, 
or infinite dimensional vector space. Let f) = (V, fj,) be a Lie algebra with Lie multiplication 
fi:V xV ^V. 

The analytic notion of contraction can be described by a continuous family of homomorphism 

U(e):V^V, £6(0,1], 17(1) = 1, 

which are nonsingular for e > and singular for e = 0. 
The new Lie bracket on V, 

He{a,b) = U(e)fi{U- 1 (e)a,U- 1 {e)b), a,b£V, e > 0, 

corresponds to a change of basis given by U(e), and leads to the Lie algebra f} e = (V, fj, £ ) 
isomorphic to f). If 

fj, (a,b) = lim/i e (a, b) 

e— >0 

exist for all a, b 6 V, we call f/ = (V, //) the contraction of f) by U(e); t)' is, in general, not 
isomorphic to f). 

For a generalized Inonii-Wigner contraction, U(e) has a simple form, namely 

JV 

y = 0y«, N>1, 

8=0 

U(e)\ v(i) = e ni t\ viih 

with 

< n < n\ < ... < n M , Ui 6 R, i = 0, ...,M. 
The necessary and sufficient condition for (|12j) to define a contraction is 

^ y « ) l/O-)) c 0y(O ) 

with 

< n, + rij, 

where all surviving structure constants are the same as for f) . The special case no 
is called a simple Inonii-Wigner contraction. 



(13) 



= and N = 1 
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We will apply a generalized Inonii-Wigner contraction to the direct sum of N copies of a 
simple Lie algebra f): 



N 



f> =f)®---©0. 

The contraction procedure will be carried out directly on the Lie-Poisson algebra associated 
with f). We denote again with {yf j^™ the coordinates on the dual of the i— th copy of g, 
satisfying {yf,yj} = Sy C® 13 yj as in formula @. 

Since we need we make the following change of basis in f): 



to 



\ ( 



\ Vn-i j 



u 2 



V N 



1 yf 

2/2 



N-l ..N-l 



a = 1, dim g, 



where the matrix defining the change of basis is a nonsingular Vandermonde matrix since we 
take Vj / Vk, 1 < j ^ k < N '. 

Note that in the new coordinates the equation (jl.3[) is satisfied if one identifies V® with the 
space generated by {yf }„™ s and n,i with i. 

We define the homomorphism U (e) by 



i = 1,...,JV- 1. 



U(e)\ v{f> = e l t\ v(i) 

We absorb the above homomorphism U(e) into the definition of the new coordinates {y~f}^=i S , 
namely hereafter we will denote 



N 



if 



The new bracket [i t is defined by 



N 

E 

fc,m=l 



k m 



{y?4} 



N 



k=l 



Cf y] +j i + j < N, 



0(e) 



i+j>N. 



After contraction we obtain the new Lie-Poisson algebra f/: 

( 



Cfy7 +j i + j<N, 







i + j>N. 



Let us give some remarks on the algebraic structure of our algebra r/. If we define the new 
vector spaces 

N-l 

^) = 0yW i = 0,...,N-l, 

j=i 

W (o) ~ ^ 
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we obtain a flag of ideals 

C W U) j>i. 

In particular 

So we have an exact sequence: 

o _> W W ^tf^Q^o 

where i is the inclusion and ir is the projection on the quotient \)'/W^; moreover we use the 
fact that tf/WM ^ y {0) ~ fl. 

We can therefore see the contraction procedure as an extension of the starting Lie algebra q 
that is known as the U N— th jet extension" of g jS]. 

In the following theorem we give a recipe to extend the contraction procedure at the level of 
the equation (J5J). 

Theorem 3.1 Let £(A) be the Lax matrix on the Lie-Poisson algebra f)' obtained from Lg(X) 
ill]) through the following procedure: 

• make the identification Aj = —ev%, 

• take a formal power series expansion in e, 

• take the limit e — > 0. 
Then £(A) satisfies 

{£(A) ® 1, 1 ® £(/*)} + [r(A - //), £(A) ® 1 + 1 ® £(/*)] = 0, (14) 
luzf/i respect to the Lie-Poisson algebra defined by 



Cfy] +j i + j<N, 
i + j>N, 



(15) 



urai/i t/ie same r-matrix given by Q). 

Proof: Setting Aj = — evi in (jllj) we easily obtain the formal expansion 

N 



L S (X; = r + P(X«) £ £ ^ (^)' ™- 



i=l j=0 

Through the identification 

N 
3=1 

we immediately get 

JV-l 

L 6 (A;e)~r+^A(A)/f(A) + 0(e), 

8=0 

where 1%(X) and ff (A) are given respectively by 

li(\)=g a pp(X a )y?. 
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Finally taking the limit e — > we obtain: 



foW, foW = r (A) 



(16) 



JV-1 



C(X) = lim Lg(X; e) = r + £ 4(A) /f (A). 

e— >0 ^— -" ' 



i=0 



Let us prove ifHjl . We set • = if i + j > N. We have: 



JV-1 



{£(\)®t,t®C(p)} = 9a/3 g-ys p(X a ) ® p{X 1 ) E ftW f](p) {yfj^i} 



i,j=0 



JV-1 



<M 5 7 5 pPH ® p(* 7 ) Cf ^ /f (A) //(/x) yl 



+3 



JV-1 



= gap g lS P (x a ) P (x-y) cf E # E //(A) fl .( M ). (17) 

i=0 j=0 

On the other hand we have: 

[r (A - /x), £(A) <g> 1 + 1 ® £(/z)] = 5 a/3 # 7?7 x 

iV-l JV-1 

x p(X«) is p{X?) f a (A - /x), (p(X^) ® 1) ^ # //(A) + (10 ppT)) ^ yV fifr) 

= gapg-Yn f a {X- p) x 

x p(x s ) ® E w + °f 7 ® E »? fiw 



8=0 



j=0 



JV-1 



JV-1 



i=0 



i=0 



JV-1 



(18) 



i=0 



-,07 



where the last expression is obtained from the previous one using the identity g lv C s 
— g^s C^ 7 and swapping the indexes 7 <->/?, a <-> 5 in the first term and 7 <-» <5 in the sec- 
ond one. 

Considering equations (fTTj) and (fTH|) we obtain 



{£(A) (8) 1, 1 ® + [r(A - M ), £(A) <g> 1 + 1 ® £(/i)] = ffa/3 s 75 p(X a ) ® p(X 7 ) x 



JV-1 



cf E 



i=0 



E ffWflM + / 5 ( A - /*)// (A) - r (a - m)//gu) 
3=0 



(19) 



Let us look in more detail at the quantity enclosed in the square brackets in equation IJ19JI . 
Using (|TB|) we have: 



Ej?ww = E^(^) to 



j=0 



i! Ua, 



1 



d 



*~3 



p(p) 



3=0 



(20) 
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and 



/ 5 (A-/x)/f(A) - / a (A -ix)f*M 
1 



rt»-ril^Ww-nA-ri^^i 



1 / d d 
i! \dA d^x 



/5(A-^)/^(A)-r(A-M)/ 5 (M) 



(21) 



Finally, from and fl^U) we obtain 

{£(A) <g> 1, 1 (8 £(m)} + [r(A - aO,£(A) ® 1 + 1 ® £(/x)] = g^g^s p(X a ) ® ppT) Cf a x 

x E # i[ (dA + d^J [/ a (A)/ 7 (M) + / 5 ( A - »)fw - n x - 
x 5 7<5 p(x a ) ® pot cf [r wr (m) + / 5 (a - m)/ /3 (a) - r (a - m)/^)] } , 

that vanishes thanks to equation Ij7]). The fact that (|15[) is a Lie-Poisson algebra is then a plain 
consequence of p4|), 

□ 



4 Extensions of su(2) Gaudin models: Lagrange chains 

In this section we will consider algebraic contractions of su(2) Gaudin models. We choose the 
matrices a 1 , a 2 , a 3 as a basis of the fundamental representation of su(2): 



ia x 1 



a 



( 



\ 



i 

1 



i(T y _ 1 
2 2 



( ^ 



V 



-1 



ia z 1 



/ 



V 



i 

-i 



where a x , a y , a z are the Pauli matrices. In this basis the r— matrix @ can be written in the 
form: 



where 



^(A) = -E/o a (A)^®a c 



(22) 



a=l 



i 

sin A 



dn 
sn 



(A) 



rational case 
trigonometric case, 
elliptic case, 



(23) 



/o 2 (A) = 



l 

sin A 

i(A) 



rational case 
trigonometric case, 
elliptic case, 



(24) 



S 



t rational case 

fo (A) = ' cot A trigonometric case, 

(25) 

f^(A) elliptic case. 

Let us consider the Lax matrix of the su(2) Gaudin magnet: 

N 3 

Lg(X) = r + E E fo(* ~ A,), (26) 

i=l a=l 

where r can be chosen as: 

• any su(2)-matrix in the rational case, 

• cr 3 in the trigonometric case, 

• the identity in the elliptic case, 

and the local variables of the model {xf }, a = 1, 2, 3, i = 1, N are the generators of the direct 
sum of TV su(2) spins with the following Lie-Poisson brackets: 



where a(3~f is the cyclic permutation of 123. 

The first natural extension of Gaudin models is obtained by contracting two copies of the 
Lie-Poisson algebra su(2) (N = 2 in the Lax matrix ()26j0 . namely 

su(2) esu(2) ~ o(4) -» e(3) 

where e(3) is the real euclidean Lie-Poisson algebra in the 3-space. The Lie-Poisson brackets of 
e(3) are given by 

{y a ,y (3 }=y\ {y a ,z?}=z\ {z Q ,^} = 0, a,/3, 7 = 1,2,3. 
After contraction we get the new Lax matrix: 

L(A) = r + ° a IfoW y a + fi (A) z a ] , (27) 

a=l 

where /q (A) is given by (123 |) .(|23 |) .l[2l) |) and /f (A) is given by Note that in the rational 

case the Lax matrix ()27() is exactly the Lax matrix of the Lagrange top ^ . This particular 
contraction procedure on the two-site su(2) Gaudin magnet has been considered in [TUj . 

Moreover we observe that the spectral invariants of the Lax matrix (|27|) are the same for 
the rational, trigonometric and elliptic cases. However this is no more true if we extend our 
Lax matrix ((2T|) to the A^-bodies case as we did for Gaudin model (see equations (|5|). H10j) . (fTT|0 . 
Namely we consider the Lax matrix associated to the Lie-Poisson algebra of the direct sum of 

copies of e(3): 

N 3 

£(A) = r + E E ffQ [/o°(A - Aj) y? + ff (A - A,) zf] . 

i=l a=l 
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4.1 The rational Lagrange chain 

The rational ©^L 
2x2 Lax matrix 



The rational (B^Li e «(3) chain, that we call rational Lagrange chain is derived from the following 



N 3 



+ 



A — A, (A -A;) 5 



(28) 



i=l a=l 

where A« € C and w € K are the parameters of the model and A £ C is the spectral parameter. 
The parameter w denotes the intensity of an external field, taken as uniform (along the chain) 
and constant (in time). 

Local variables of the model are generators of the direct sum of N e(3) tops, yf,zf, i = 
1, N, a = 1, 2, 3 with the following Lie-Poisson brackets: 

{yf, yf } = 6a y], [yf, zf} = k 3 z], [zf, z?) = 0, (29) 

where a(3"f is the cyclic permutation of 123. 

The Lie-Poisson brackets (|29|) have 2N Casimir functions: 

a=l a=l 

As we have shown in the previous section the Lax matrix (|28|) satisfies the linear r-matrix 
algebra 

{£(A) ® 1, 1 ® £(fi)} + [r(A - ft), £(A) ® 1 + 1 ® £(//)] = 0, (30) 
where 1 is the 2x2 identity matrix and the r-matrix is given by 

3 



r(\) = -^a a ®a«. (31) 

a=l 

The spectral curve T, 

r : det(£(A) - /xl) = 0, 
is an hyperelliptic curve of genus g = 2N — 1, given by 

D Q ^(1) ^(2) 

-ft* <->?: 



-A-Ai (A -A,) 2 (A -A;) 3 (A-A,) 4! 
with the Hamiltonians Ri and Si given by: 



* " A"^ + 



, 44 + gfj + yf4 - vl4 ~ ^ 2 4 ~ 2 44 + 44 + 44 

(A, - A,) 2 (A, - A fc ) 3 

s. = ^ [(4) 2 + (y 2 ) 2 + (y?) 2 ] W + 

, y f 44 + gig? + g|4 | 44 + 44 + 44 

h\ x i~ x k (Ai-A fc ) 2 
These are integrals of motion of the rational Lagrange chain, which are Poisson commuting: 
{Ri,Rj} = {S h Sj} = {R t , Sj} = 0, i,j = 1, N. 
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Notice that there is one linear integral given by 

N N 



i=i 



i=i 



Let us remark that if N = 1 the Lagrange chain degenerates into the well-known e(3) Lagrange 
top, whose Hamiltonians are given by y 3 and (y 1 ) 2 + (y 2 ) 2 + (y 3 ) 2 + 2wz 3 . 
We can bring the curve T into the canonical form by the scaling 



N 



p\ — *p = n 



i=l 



The equation of the spectral curve becomes 



w 



N p <j r (i) 

2 + E^+ - " + 



+ 



(2) 



-A- A, (A — Aj) 2 (A - A;) 3 (A -A;) 



v 



i=i 



10 A + SiA + S2A + ... + S47V, 



(32) 



where the coefficients s 3 -, j = l,...,4A r are linear combinations of the Hamiltonians and of the 
Casimir functions. 



4.2 The trigonometric Lagrange chain 



The trigonometric Cj(3) chain, that we call trigonometric Lagrange chain, is derived from 

the following 2x2 Lax matrix 



N 3 



£(A) = wa 3 + £ E ° a [ /o°(A - A,) yf + /f (A - A,) zf] 



(33) 



i=l o=l 



where 



1 



sin A 
cot A 



/o 2 (A) 
/i 2 (A) 



1 



sin A 
cot A 



/ 3 (A) = cotA, 
/f(A) = cot 2 A + 1. 



sin A ' J 1 v ' sin A 

The Lax matrix (|33[) satisfies the linear r-matrix algebra (j3()[) with r-matrix given by (|22|24|) 

The spectral curve V, 

r : det(£(A)-/xl) = 0, 

is given by 



-A* 



Q + ^cot(A - Xi)Ri + cot 2 (A - A*) 5, + cot 3 (A - Ai) cf } + cot 4 (A - A 4 ) cf \ 



,(2) 



i=l 



where 



N 



Q 



w 



+ ^2[2wzf + (y}) 2 + (y 
1 



2\21 



i=l 
AT JV 

i=l fc^i 

2 

+ 



sin (Aj - A fc ) 



• . , x x [yl4 + yhf + cos(\ l - x k )ylzf] + 

sin^A, — Ak) 

[cos(Ai - A fc ) (4^ + ) + *ik 3 ] - 



11 



k^i 



N 

,l„.l i „.2„,2 



sin(Ai - A fc ) 



VkVi + VkDi + cos ( A < - x k)yM- 



sin 2 (A; - A fc ) 



sin 2 (Ai - A A 



Z k 2 



+ 



l 



sin 2 (Aj - A fc ) 
l 



[cos(Ai - X k )(yl4 + yf4 - v\ z \ - y\ z h + y 4 3 4 - yhf] 



(yl) 2 + ivif + (yf) 2 + + ft)* + {zff + (z?y + 

N 



sin(Aj - A fc ) 



x z\ + y 2 ,^ + cos(Ai - A fc )y^ 3 ] 



sin 2 (Ai 



— [cos(Ai - \k){z\z\ + zlzf ) + zlzf] \ . 



Note that the coefficient Q is functional dependent on the Si, Ri and C f 

2 /AT /ON \ 2 



(2). 



N ( N r \ ( N 
i=i \i=i / \j=i / 



As in the rational case we identify in the coefficients Ri, Si the Poisson commuting integrals 
of the trigonometric Lagrange chain: 



{Ri,R 3 } = {Si,Sj} = {Ri,Sj} = 0, i,j = l,...,N. 



5 Concluding remarks 

In this paper we have performed a Inonii-Wigner contraction on Gaudin models, showing that 
the integrability property is preserved by this algebraic procedure. Other interesting structures 
can be inherited from the Gaudin models as well. For example in JOJ it is shown how to obtain 
Backlund transformations for the Lagrange top from those of the su(2) Gaudin model. It turns 
out that it is possible to generalize this approach to the construction of Backlund transformations 
for Lagrange chains. This work is in progress and we plan to report the results in a forthcoming 
paper. 

Finally we want to remark that the contraction procedure described in this paper can be 
applied to quantum Gaudin models as well, in a straightforward way. 
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Appendix: The general rational Lax matrix 

In the proof of the theorem 3.1 we have seen that our Lax matrix £(A) satisfies the r-matrix 
structure 

{£(A) 81,10 C(jjl)} + [r(A - /j), £(A) ® 1 + 1 ® C(p)] = 
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provided that the functions /"(A) satisfy the following system of functional equations (see equa- 
tion (dU): 



dim g 

13,5=1 



£ » + / 5 (A - /x)/f (A) - r (A - M )//0*) 

i=o 



(34) 



for all a, /?, 7 and i = 0, . . . , N — 1. 

In theorem 3.1 we have shown that the system (|34|) is satisfied with the choice IjlfiJ): however 
this is not the general analytic solution to the system (|34[). In this appendix we give a conjecture 
about the general analytic solution of (|34jl in the rational case, where we have 

/o a (A) 



A - c ' 



Co being an arbitrary constant; /^(A) obviously does not depend on a. 
The system of functional equations (|34j) reduces to the following one: 



i=0 



7fc(A) - 

[A — \ 



k = 0,. 



(35) 



Our solution is based on the following algebraic conjecture: 



Conjecture 1 Let k G N, c\, . . . , c& arbitrary constants and let us denote with {Jt} the set of 
vectors satisfying the Diophantine equation: 

{J»} = {JGN i : J 1 +2J 2 + ---+iJ i = i}, 

then the following polynomial identity holds: 



EIE^-^ 17 " E V- 



c/ 1 c Jfc 



t-1 <-l. 



i=0 {-4-*} 



Ji! " Jfc-i! 
where we denoted with |Jj| £/ie 1—norm of J{: 

i 

\Ji\ = Y. Ji - 
1=1 

Using the program for symbolic manipulation Maple we tested this conjecture for N < 20. 

Using the explicit form for the function /o(A) and our conjecture it is not difficult to construct 
the general analytic solution for the system (|35[). reading: 



Ji 



/o(A) 



l,...,iV. 



(36) 



13 



In fact, inserting (|3fi|) into (|33|) we get: 



k 



£/i(A)/ fc -i(M) 



i=0 





\dX dfi J \ [j, — X 



) 



H — A 



The fact that solutions (|36|) provide the general analytic solution to (|35|) can be argued taking 
the limit A — ► fx in the functional equations (|35|) . yielding the system of ordinary differential 
equations: 



As the system (|37(l is triangular, for any given i the functions fi(X) can be found solving a 
system of i + 1 ordinary differential equations and therefore depend at most upon i + 1 arbitrary 
parameters. This is exactly the number of arbitrary parameters entering our solutions (|36|). 

Using our conjecture it is straightforward to extend the solutions (|36[) to the trigonometric 
and elliptic case. It is enough to replace fi(X) with /"(A) and /o(A) with the corresponding 





(37) 



/ Q (A) given by (|23|24|25|): 




i = l,...,N. 
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